Near- and mid-infrared photometry of high-redshift SCR sources 
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ABSTRACT 

Using the Spitzer Space Telescope, we have obtained 3.6-24 /im photometry of 38 radio galax- 
ies and 24 quasars from the 3CR catalog at redshift 1 < z < 2.5. This 178 MHz-selected sample 
is unbiased with respect to orientation and therefore suited to study orientation-dependent ef- 
fects in the most powerful active galactic nuclei (AGN). Quasar and radio galaxy subsamples 
matched in isotropic radio luminosity are compared. The quasars all have similar spectral energy 
distributions (SEDs), nearly constant in i/Fi, through the rest 1.6-10 ^m range, consistent with 
a centrally heated dust distribution which outshines the host galaxy contribution. The radio 
galaxy SEDs show larger dispersion, but the mean radio galaxy SED declines from rest 1.6 to 
3 iim and then rises from 3 to 8 /xm. The radio galaxies are on average a factor 3-10 less lumi- 
nous in this spectral range than the quasars. These characteristics are consistent with composite 
emission from a heavily reddened AGN plus starlight from the host galaxy. The mid-infrared 
colors and radio to mid-infrared spectral slopes of individual galaxies are also consistent with 
this picture. Individual galaxies show different amounts of extinction and host galaxy starlight, 
consistent with the orientation-dependent unified scheme. 

Subject headings: Galaxies: active - quasars: general - Infrared: galaxies 



1. Introduction 

When exploring the general evolution of galaxy 
populations across cosmic times, a particular chal- 
lenge is to distinguish between black hole and star- 
forming activity. Star formation and obscuring 
dust go hand in hand, and black-hole-driven ac- 
tive galactic nuclei (AGN) are also surrounded by 
dust mainly distributed in a disk/torus-like ge- 
ometry (Antonucci 1993). There is evidence that 
AGN mainly power the near- and mid-IR emis- 
sion (NIR, ^2/im; MIR, ^10 iim) from hot nuclear 
dust, while starbursts contribute mainly to the far- 
infrared (FIR, ~60 /im) luminosity (e.g., Rowan- 
Robinson 1995, Vernet et al. 2001, Schweitzer et 
al. 2006). Using the MIR/FIR luminosity ratio as 
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an indicator for the relative AGN and starburst 
contributions, numerous studies have found an in- 
crease of AGN/starburst activity with total lumi- 
nosity and redshift, but the validity of this trend is 
still under discussion because of selection effects on 
the various samples. More seriously, an unfavor- 
able AGN orientation could cause MIR obscura- 
tion (e.g. Pier & Krolik 1993), leading to a funda- 
mental observational degeneracy: a low MIR/FIR 
luminosity ratio can be due to either a high star- 
forming contribution or to an AGN in which the 
hot dust is obscured. The spectral energy dis- 
tribution (SED) of an obscured AGN may thus 
mimic that of a starburst-powered source. While 
this degeneracy has now been widely examined at 
low/intermediate luminosity and redshift {z < 1), 
it has still to be explored for the most luminous 
sources at high redshift (z > 1). In order to as- 
sess galaxy and AGN evolution in the universe, we 
therefore need to understand this AGN/starburst 
degeneracy for a population of luminous high- 
redshift sources. A crucial step towards this is to 
study the orientation dependence of the NIR and 
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MIR emission of high-redshift AGN. 

Orientation-dependent effects can only be 
tested and quantified with AGN samples hav- 
ing type 1 (unobscured) and type 2 (obscured) 
subsamples matched in isotropic emission. The 
clean AGN tracers — optical, [O III] A5007A, 
NIR, and X-ray (<10 keV) — all fail to fulfill 
this requirement. The [O II] A3727A emission, 
while isotropic (Hes et al. 1993), is probably dom- 
inated by extended starbursts and shocks (Best 
et al. 2000) rather than by the AGN. Therefore, 
the only feasible way is low-frequency (meter- 
wavelength) radio selection because the integrated 
emission from the radio lobes is optically thin and 
essentially isotropic. This makes radio-loud AGN 
particularly attractive for studying orientation- 
dependent properties at other wavelengths and, 
after sorting out the influence of radio jets/lobes 
on the emission, for generalizing conclusions about 
orientation-dependent effects to the much larger 
population of radio-quiet AGN. 

The brightest low-frequency-selected AGN 
sample is the SCR compilation (Spinrad et al. 
1985). The powerful double- lobed radio galaxies 
(henceforth simply called radio galaxies) are sup- 
posed to be misaligned quasars (Barthel 1989). 
Based on IRAS coadded scans and a few indi- 
vidual detections, Heckman et al. (1992, 1994) 
already noted an average MIR/FIR difference be- 
tween 3CR quasars and radio galaxies. More com- 
prehensive MIR and FIR spectrophotometry from 
ISO is in hand (as compiled by Siebenmorgen et 
al. 2004 and by Haas et al. 2004) as well as from 
Spitzer (e.g., Shi et al. 2005, Haas et al. 2005, 
Ogle ct al. 2006, Cleary et al. 2007), providing 
a basis to study the z < 1 SCR objects. These 
sources are, however, a factor of five less radio- 
luminous on average than the most powerful radio 
sources seen at higher rcdshift, and the lower indi- 
cated accretion power may reflect diflerent source 
physics. The higher-luminosity population can 
be sampled by the SCR sources at 1 < z < 2.5, 
which have radio luminosities similar to those of 
the most powerful radio sources at even higher 
redshift (2.5 < z < 6). However, with the excep- 
tion of a few targets (Siebenmorgen et al. 2004, 
Seymour et al. 2007), the high- 2; SCR sample has 
not been well observed in the rest frame NIR and 
MIR. 

This paper is based on Spitzer observations of 



62 of the 64 high-z SCR sources. It focuses on the 
hot nuclear dust emission and its obscuration in 
the most luminous type-2 AGN. We use a ACDM 
cosmology with Hq = 71 km/s/Mpc, Om = 0.27 
and = 0.73. 

2. Observations and Data 

With the Spitzer Space Telescope (Werner et 
al. 2004), we are observing the entire sample of 64 
high-z SCR sources using the instruments IRAC 
(3.6-8.0 /xm, Fazio et al. 2004), the IRS blue peak- 
up array (16 /um, Houck et al. 2004) and MIPS 
(24 /zm, Rieke et al. 2004). Most observations 
are performed in our guaranteed time program 
(PID 40072; PI G. Fazio) with on-source exposure 
times 4x30s (each IRAC band), 4xl4s (IRS), and 
10 X 10s (MIPS). A few sources have been observed 
in other programs, and we use the published pho- 
tometry if available (e.g., PID 3329; PI D. Stern, 
Seymour et al. 2007). 

For IRAC, we used the basic calibrated data 
products (BCD, version S16) and coadded them to 
0'.'869 pixels using the latest version of IRACProc 
(Schuster et al. 2006) . This optimally handles the 
slightly undersampled IRAC PSF in order to as- 
sure accurate point-source photometry. For IRS, 
we used the post-BCD pipeline product, version 
S16. For MIPS, we used custom routines to mod- 
ify the version S17 BCD files to remove instrumen- 
tal artifacts (e.g., residual images) before shift- 
ing and coadding to create the final mosaics. All 
sources are well seen on the images in all filters. 
The sources were extracted and matched using the 
SExtractor tool (Bertin & Arnouts 1996). We used 
sufficiently large apertures so that aperture cor- 
rections are small (<5%). The photometric errors 
are typically smaller than 10% but increase for 
faint sources; exceptions are 3C 225A and 3C 294, 
where nearby bright stars make the photometry 
uncertain in the shorter IRAC bands. 

As of 2008 April, 24 quasars and 38 radio 
galaxies have been observed, covering the com- 
plete high-z SCR sample with the exception of 
the quasar 3C245 and the radio galax;y 3C325. 
All 62 sources have IRAC measurements and are 
observed in at least one of the 16 and 24 fim bands 
(54 sources at 16 /um and 60 sources at 24 /zm). 

For the analysis it is desirable to compare rest 
frame SEDs with the same wavelength sampling. 
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Fig. 1. — Infrared versus radio luminosity of the 3CR sample at 2 > 1 prior to normalization, 'x' symbols 
denote quasars; circles and squares denote radio galaxies. Superposed crosses indicate radio galaxies that 
show evidence of silicate absorption ( §3.ip . The vertical long-dashed lines mark the range of our luminosity- 
matched quasar and radio galaxy subsamples. The dotted lines indicate Ls^m/ LnsMUz ratios of 1, 10, and 
100. The radio galaxies are grouped into several SED classes in Fig.[3]and ij31l The color-coding and symbols 
are: green circle (A), red circle (B), red square (C), blue square (D), blue circle (E). The two low-excitation 
radio galaxies 3C68.1 and 3C 469.1 are labeled with their 3C numbers, as are sources outside the luminosity 
range of our analysis. 



Depending on the redshift of our sources (1 < z < 
2.5) the observations sample different rest wave- 
lengths between 1.6 and 10 /im. Before resampling 
and interpolating the SEDs, we checked that spec- 
tral features do not affect the interpolation. In 
principle, prominent spectral features in this wave- 
length range could be PAH emission bands around 
7.7 fiva and the 9.7 /xm silicate absorption feature. 
IRS spectra of several \ow-z and even a few high- 
z PR II radio sources are available. PAH features 
are weak and usually undetected, and the continua 
are generally smooth. Strong silicate absorption 
is, however, present in some objects (Haas et al. 
2005, Ogle et al. 2006, Cleary et al. 2007, Sey- 
mour et al. 2008). Our broadband SEDs therefore 



represent the smooth continua for A < 8 /im but 
are uncertain at rest wavelengths near 10 /im. In 
practice, the de-redshifted SEDs were interpolated 
in log-log space at 12 wavelengths between rest 1.6 
and 10 /im to produce the figures shown. 

The quasar and radio galaxy samples match 
reasonably well in redshift and rest-frame 178 MHz 
flux density. Rest-frame 178 MHz radio flux den- 
sities were derived from data listed in the NASA 
Extragalactic Database, NED. The quasars have 
mean redshift (z) = 1.44 ± 0.31 and mean flux 
density (Sns) — 27.8 ± 15.1 Jy, while the val- 
ues for the radio galaxies are (z) — 1.42 ± 0.31, 
(Sns) = 22.2 ± 6.2 Jy. Thus over the whole 
sample, quasars are about 30% more luminous 
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Fig. 2. — Rest frame quasar and radio galaxy SEDs normalized by rest 178 MHz luminosity {{Lns) ~ 
5.4x10^'* erg/s). Symbols connected with thick blue and red lines show the mean SEDs for quasars and 
radio galaxies, respectively. The thin dotted blue and red lines indicate the dispersion (upper and lower 
quartiles) around the mean SEDs; the mean ratio of upper/lower quartiles are 2.3 (quasars) and 3.4 (radio 
galaxies). The radio galaxy SED can be explained by the sum (black long-dashed line) of a reddened quasar 
(blue long-dashed line) and starlight from the host galaxy (thin black solid line) . The long-dashed lines have 
been shifted slightly to make them visible in the plot. The difference between radio galaxies and reddened 
quasars at 10 /im may be due to the silicate absorption feature which may escape detection in our broad 
band photometry. 



than radio galaxies as shown in Figure [TJ In or- 
der to improve the luminosity match, we have 
excluded the sources at the low and high ends 
of the distribution, Lns < 1-8 x 10^^ erg/s and 
Lns > 1.5 X 10*^ erg/s, respectively. We have 
also excluded the quasar 3C418 because of its 
flat radio spectrum (low-frequency spectral index 
ans « 0), while all other sources have steep ra- 
dio spectra (—1.1 < ans ^ —0.6). The resulting 
mean radio luminosities of the sample galaxies are 
(Lirs) = (5.35 ±2.53) x 10^"* erg/s for quasars and 
(5.55 ± 2.34) X 10'''* for radio galaxies. 

While the quasar and radio galaxy distribu- 
tions match very well in L^s, a proper analy- 



sis of orientation-dependent effects requires also 
that the individual SEDs are normalized by the 
radio luminosity, which serves as a tracer for the 
intrinsic AGN strength. Therefore, we have nor- 
malized each SED to the sample mean 178 MHz 
luminosity; after normalization each object has 
Lns = 5.4 X 10^"* erg/s. Because of the good Lns 
match of the samples, it turned out that the net 
effect of this normalization on the results is small. 
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3. Results and Discussion 

3.1. Radio galaxies as obscured quasars 

The NIR-MIR SEDs of quasars are all very sim- 
ilar in shape, as shown in Figure [2l The SEDs 
can be described by a single power law cx v^^, 
consistent with previous results for lower-redshift 
objects (e.g., Elvis et al. 1994). The dispersion 
of the SEDs is essentially caused by differing ra- 
tios of MIR to radio luminosity. Some quasars ex- 
hibit small (10-20%) bumps around 5 um explain- 
able by distinct hot dust componentso The power 
law shape of the quasar SED can naturally be ex- 
plained by the superposition of centrally-heated 
dust components with a radial temperature gra- 
dient (1500 K > T > 300 K) as has been found 
also in lower luminosity type-1 AGN (e.g. Ward et 
al. 1987, Barvainis 1987; see also Rowan-Robinson 
1980). Any contribution of the quasar host galax- 
ies to the NIR-MIR SEDs appears to be outshone 
(factor >5-10) by the AGN dust emission. 

In contrast to quasars, radio galaxies display a 
diversity of SED shapes leading to a 50% larger 
dispersion around their mean SED (Fig. De- 
spite the dispersion, nearly all radio galaxy SEDs 
show a decline from rest 1.6 /um to 3 /xm and a 
rise from 3/im to 8 /zm (L^ oc In addi- 

tion, the average radio galaxy SED is fainter by 
a factor of three at 8 fim and a factor of eight 
at 2/Ltm relative to the quasar SED. Unlike the 
quasars, hot {T > 750 K) dust emission is not 
seen in the radio galaxy SEDs. Its absence can 
be explained by absorption (screen Ay ~ 50 jl of 
the central dust emission. The short wavelength 
(A < 3 ^m) component can then be explained by 
emission from stars in the host galaxy. Extrap- 
olation of the mean 3-8 fim SED slopes towards 
longer wavelengths suggests that the radio galaxy 
and quasar SEDs meet each other at about 25- 



^ Despite the similarity of the infrared SEDs, the quasar pop- 
ulation is not homogeneous at optical wavelengths: there 
are quasars like 3C 186 with blue optical SED and 3C68.1 
with red optical SED, as listed in NED. In the orientation- 
based unified scheme, 3C68.1 could be borderline so that 
the broad lines are detected, but most of the UV-optical 
continuum is absorbed. 

^ The reddening curve used is a compromise be- 
tween the latest results for Milky Way reddening 
and earlier data (summarized by Indebetouw et 
al. 2005): Ay / Afl- / As^.^ / As^m / Ag^m / Aio^im = 
1 / 0.184 / 0.070 / 0.037 / 0.028 / 0.040. 



40 /zm, and beyond these wavelengths extinction 
may be no longer relevant. 

As noted above, the quasar NIR-MIR SED 
shapes are extremely homogeneous. This is re- 
flected in the narrow range of the quasars' NIR 
and MIR colors. The color-color diagram shown in 
Figure[3]illustrates the differing SED types. In this 
diagram, quasars populate a distinct locus ("E"), 
while radio galaxies show wider dispersion as men- 
tioned above. According to their location in the 
color-color diagram, we have grouped the radio 
galaxies into five classes described below. Their 
SED shapes are illustrated in Figured! 

A) Four sources at the high end of the 3/im/1.6/im 
ratio (above the dotted and dashed boxes 
in Fig. [3]): basically, they have quasar- 
like SEDs, but the hottest dust emis- 
sion at about 1-2 /im appears to be ab- 
sorbed (screen Ay « 5) leading to a redder 
3/Ltm/1.6 /xm color compared to quasars. 

B) The bulk of radio galaxies (20 sources) have 
declining 1-3 /xm SEDs with a steep 3-8 ^m 
rise. Their colors can be explained by a 
heavily reddened AGN plus an added com- 
ponent of starlight contributing at 1.6 fim. If 
this explanation is correct, the direction of 
the extinction vector Ay becomes meaning- 
less because host galaxy starlight will not 
be affected by extinction near the nucleus. 
Instead, vertical position in the plot is de- 
termined by the relative contributions of 
starlight and AGN light, while horizontal po- 
sition measures the amount of extinction (to 
the extent the underlying AGN SEDs are the 
same). As noted above. Ay ~ 50 mag is re- 
quired to explain the colors. 

C) Three sources at the low end of the 3 /im/1.6 fim 
ratio (below the long dashed line in Fig. [3]): 
Their SEDs show a very strong host galaxy 
contribution at 1.6 /zm, and starlight exceeds 
the dust luminosity even at wavelengths as 
long as 3.5 /im. In principle, class C is simi- 
lar to class B but with stronger host galaxy 
contribution. 

D) Three sources immediately below the dot- 
ted box in Fig.m Their SEDs can be ex- 
plained by a slightly reddened AGN (sim- 
ilar to class A) plus an added component 



5 











CO 

^ 4 


■ - 


A ° 

?...d267'^ 


\ J, 469 1 

\ 


E 




• > X 


X X 1- 


CD 




B * .• 

• • , -^68.2 ■ 

L • '< ^ 


D ~: 


E 


_ 






F3 




■ ■ 


X quasars 
° o galaxies 
+ Si absorption 






0.3 0.5 


0.8 



F 5 |im / F 8 |Lim (rest) 



Fig. 3. — NIR/MIR color-color diagram. The radio galaxies are grouped into five classes labeled A-E as 
explained in §3.11 Symbol color-coding is the same as in Figure [T] Radio galaxies with a photometric 
signature for silicate absorption are additionally marked with an underlying plus. The two sources 3C267 
and 3C 469.1 with spectroscopically-detected silicate absorption are labeled, as are the two low-excitation 
radio galaxies 3C68.1 and 3C 469.1. The error bar in the upper left corner represents a color rms of 15%. 
The Ay arrow indicates screen extinction with the reddening law given in i )3.1l 



of starlight contributing significantly at 1.6- 
3 ^m. 

E) Three sources with quasar-like SED colors 
(inside the dotted box in Fig.g]): Their SEDs 
overlap with the low luminosity end of the 
quasar SEDs. In the orientation-based uni- 
fied scheme, these sources could be border- 
line so that most of the dust torus is visible 
but the broad line region and the UV-optical 
continuum are obscured. 

While the rest 8-10 /im range is poorly sam- 
pled, eight galaxies show declines in this range 
that could be caused by sihcate absorption. (The 
MIPS-24 filter, 50% transmission at 20.8-29.3 ^m, 
requires z < 1.8 for the silicate feature to fall 
into its range.) One of these sources (3C 469.1, 



z = 1.336) has an IRS spectrum available. It 
shows a broad silicate absorption with optical 
depth rg.7 « 0.55 corresponding to Ay ~ 10, con- 
sistent with its position in Fig. [3] This supports 
the view that the SED declines in the other ra- 
dio galaxies are due to silicate absorption, too. 
The photometric silicate absorption sources show 
a wide range of colors (Fig. [3]), but only one galaxy 
(3C 469.1) is on the blue (right) side, where low- 
extinction sources reside. The important conclu- 
sion is that the silicate feature requires consider- 
able extinction to be present in at least some of 
the radio galaxies, and this is largely independent 
of the SED classi 

^ The photometric silicate absorption sources are 3C68.2, 
3C222, 3C249, 3C250, 3C266, 3C 305.1, 3C324, and 
3C 469.1. These galaxies lie in the redshift range 1.08 < 
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Fig. 4. — Mean SED of each radio galaxy class as identified in The SEDs have been normalized to 

the mean f 78 MHz luminosity. The mean quasar SED is also shown for comparison. The dispersion around 
each SED (measured as mean ratio of upper/lower quartile) is: 2.3 (quasars), 2.0 (A), 3.2 (B), 2.4 (C), 1.7 
(D), 1.5 (E), and 4.8 (sihcate absorption). 



If radio galaxies are misaligned quasars, as pro- 
posed in the unified scheme, reddening of individ- 
ual galaxies should be correlated with their ex- 
tinction. Figure [5] shows that this is indeed the 
case. Quasars populate a distinct region of this di- 



2 < 1.83, suggesting that in this range the broad band 
16/im/24/im filter combination is able to register silicate 
absorption, if strong enough. For comparison, this redshift 
range contains 20 more radio galaxies with 16 and 24 fim 
photometry available but without silicate absorption sig- 
natures in their broadband SEDs. Four of these sources 
(3C 13, 3C 266, 3C 267, 3C 356) have IRS spectra available, 
but significant silicate absorption is detected only in one 
of them (3C267, z = 1.14, rg.r ~ 0.2). At low redshift, 
the detection of silicate absorption appears not to be di- 
rectly correlated with other absorption signatures, perhaps 
because of complex geometry and/or varying silicate dust 
abundance (e.g., Haas et al. 2005, Ogle et al. 2006, Cleary 
et al. 2007). A detailed analysis of the high-2 3CR spectra 
and the photometric detectability of silicate features will 
be presented elsewhere (Leipski et al. in prep.). 



agram characterized by high MIR/radio and blue 
NIR-MIR colors. Most radio galaxies spread to- 
wards fainter MIR/radio and redder NIR/MIR. 
Under the reasonable assumption that the emis- 
sion at 5-8 /im is not affected by the host galaxy, 
de-reddening along the direction of the extinction 
vector can place each radio galaxy inside the re- 
gion populated by quasars. Thus individual ra- 
dio galaxies can be explained as reddened quasars, 
consistent with the orientation-dependent unified 
scheme. 

The typical amount of radio galaxy reddening, 
Ay ~ 50 for an obscuring screen (Fig. [5]), cor- 
responds to a hydrogen column density Nh ~ 
9 X 10^^ cm^^ (for Ay — 5.6 x 10^^^ mag/cm^^ 
— Seward et al. 1999). This is close to but be- 
low the Compton-thick limit {Nh = 10^** cm~^). 
Screen extinction is a simplification, and one may 
expect a more complex geometry. If emitting dust 
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Fig. 5. — Infrared-radio color-color diagram. The dotted line marks the region occupied by quasars. The 
color-coding and symbols of radio galaxies correspond to those in Figures [T] and [3l The two sources 3C267 
and 3C 469.1 with spectroscopically-detected silicate absorption are labeled as well as the two low-excitation 
radio galaxies 3C68.1 and 3C 469.1. The error bar in the upper left corner represents an rms of 15%. The 
Ay arrow indicates screen extinction with the reddening law given in tJ3.1l 



particles are spatially mixed with the absorbing 
ones, the amount of dust has to be higher for 
the same observed reddening, typically by a factor 
3-50 Thus there could very well be enough gas 
present to render the AGN Compton-thick. 

There is, unfortunately, no independent mea- 
surement of reddening for individual galaxies, nor 
is it certain that a Galactic reddening curve ap- 
plies to AGN. Thus it is still an open question 
whether after de-reddening there will remain a dif- 
ference in the 8 fim/178 MHz ratio between radio 
galaxies and quasars. If such a difference remains, 
with quasars having a higher 8 /im/178 MHz ratio 
than radio galaxies, then either our screen extinc- 



* The transmission factors are exp{—T\) and (1 — 
eKp{—Tx)) /tx for the screen and the mixed case, respec- 
tively, with Tx = 0.916 X Ax (Disney et al. 1989). 



tion premise is too simple or the MIR luminosity 
of quasars is enhanced by an additional — po- 
tentially non-thermal — contribution. Our Chan- 
dra X-ray observations of a subset of the sample 
will provide independent estimates of the extinc- 
tion towards the nuclei (Wilkes et al., in prep.). 

To summarize, while quasars exhibit a uniform 
SED shape which can be explained by a centrally 
heated dust distribution, radio galaxies show a di- 
versity of SED shapes. In all cases, however, the 
radio galaxy SEDs are consistent with being in- 
trinsically a quasar modified by absorption of the 
dust emission and addition of some amount of host 
galaxy starlight. 



8 



3.2. Evolutionary effects and non-thermal 
contributions 

Studying powerful 3CR sources at z < 1, 
Ogle et al. (2006) found evidence for a popu- 
lation of accretion-inefficient radio galaxies, in 
which the jet/lobe may be powered by extrac- 
tion of rotational black hole energy. These 
sources, mainly optically-classified low-excitation 
radio galaxies (LERGs), have a 15 /im luminos- 
ity below 8 X 10^'^ erg/s and a luminosity ratio 
ii5/^m / Li78MHz < 10. In Contrast, with the reaon- 
able assumption that Lg^m ^ ii5/^m, all our 
z > 1 sources have observed MIR luminosity 
Lsiim > 5 X 10''^ erg/s, which is expected to be 
even higher after de-reddening. Also, the two 
LERGs (3C68.2, 3C 469.1) in our sample show a 
high luminosity ratio Lg^jm / Litsmhz > 10 compa- 
rable to that of quasars (Fig.[T]). From this, our 
data do not support the existence of an accretion- 
inefficient population among the powerful SCR 
sources at z > 1. A possible explanation for the 
deficit of optical high-excitation line luminosity 
(for instance [O III] A5007A) in our two LERGs 
may be extinction of the narrow-line region. On 
the other hand, some of our radio galaxies with 
very strong host contribution (plotted as squares 
in Fig.E]) are expected after de-reddening to lie 
at the low end of the Lg^m/^iTSMHz distribution. 
Hence compared with the strength of both the host 
and the radio lobes, these galaxies are relatively 
weak in the MIR and may represent a population 
at the beginning of a different evolutionary state. 

Some authors have attributed the excess emis- 
sion of quasars compared to radio galaxies to non- 
thermal emission from synchrotron jets. For ex- 
ample, Cleary et al. (2007) fitted the SEDs and 
spectra of SCR sources at 0.5 < z < 1 with a com- 
bination of a spherically symmetric dust model 
and a jet-|-lobe synchrotron component. They at- 
tributed half of the factor of four excess in the 
15 /im luminosities of steep-spectrum quasars rela- 
tive to radio galaxies to a non-thermal component. 
If such a non-thermal component were present in 
our SCR sources at z > 1, it would show up in 
Fig.[5]as an offset by about a factor of two between 
dereddened radio galaxies and quasars. This con- 
clusion is, however, dependent on both the red- 
dening law and on the radiative transfer and thus 
the geometry of the emitting region. In order to 
draw definite conclusions about any MIR luminos- 



ity excess, detailed radiative transfer modelling is 
required (Heymann et al., in prep.). Spherically 
symmetric models are wholly inadequate for this 
purpose. In an inclined disk-like system, some 
fraction of the MIR emission is likely to have very 
little obscuration while the bulk of the MIR emis- 
sion is heavily obscured, and no spherical model 
can properly account for this geometry. All we can 
say at the moment is that our data appear con- 
sistent with a simple thermal interpretation and 
show no evidence for a non-thermal component. 

4. Conclusions 

The SCR sample at 1 < z < 2.5 represents 
the most luminous steep-spectrum quasars (type 1 
AGN) and powerful double-lobed radio galaxies 
(type 2 AGN). This sample is nearly unbiased by 
orientation. We have defined subsamples of 19 
quasars and SS radio galaxies matched in isotropic 
rest 178 MHz luminosity and have obtained Spitzer 
S.6-24 fim photometry. The main results are: 

1) Quasars all have similar energy distributions 
in the rest frame 1.6-10 /im range, and their 
ratio of MIR to radio luminosity is also 
nearly constant. This is consistent with re- 
sults seen previously in lower-redshift sam- 
ples. 

2) The rest frame 1.6-10/im SEDs of ra- 
dio galaxies can be explained as reddened 
quasars, consistent with orientation-dependent 
unification. Various amounts of extinction 
of the AGN emission combined with addi- 
tion of host galaxy starlight can explain the 
diversity of radio galaxy SEDs. 

S) If the extinction is sufficiently large, there 
is no need to invoke a beamed synchrotron 
contribution to explain the MIR luminosity 
difference between quasars and radio galax- 
ies. The actual amount of extinction has to 
be derived from additional observations. 

4) The above results hold also for splitting our 
sample in redshift and luminosity; within 
our sample we do not find any trends with 
redshift or luminosity. 

5) At rest frame 8 /im, quasars are S times more 
luminous than radio galaxies. If this dif- 
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ference applies also to high-redshift, radio- 
quiet AGN, then MIR (24 /xm) surveys are 
strongly biased in favour of type-1 and 
against type-2 AGN. This will make it very 
difficult to resolve the AGN/star burst de- 
generacy with only broadband SEDs, and 
spectral line diagnosis will be required. 

While our ncar-mid-IR SEDs provide a funda- 
mental set of high luminosity AGN templates, we 
expect to derive more stringent conclusions from 
proper two-dimensional radiative transfer mod- 
elling in combination with Spitzer MIR spectra, 
Chandra X-ray observations, and Herschel far- 
IR/sub-mm data. 
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